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2-color QCD at finite temperature T and chemical potential [i is revisited in the strong 
coupling limit on the lattice. The phase structure in the space of T, n and the quark mass 
m is elucidated with the use of the mean field approximation and the large dimensional 
expansion. We put special emphasis on the interplay among the chiral condensate, the 
diquark condensate and the quark density in the T-/j,-m space. Qualitative comparison is 
made between our results and those of recent Monte-Carlo simulations in 2-color QCD. 

§1. Introduction 

Physics of matter under high baryon density is one of the most challenging 
problems in quantum chromodynamics (QCD) both from technical and physical 
point of view. Unfortunately, analysing the finite baryon chemical potential region 
of 3-color QCD from the first principle lattice simulations is retarded due to the 
complex fermion determinant. However, the situation is different for 2-color QCD in 
which the fermion determinant can be made real and the Monte-Carlo simulations 
are attainable. Therefore, 2-color QCD provides an unique opportunity to compare 
various ideas at finite chemical potential with the results from lattice simulations. 

One of the major differences between 2-color QCD and 3-color QCD is that the 
color-singlet baryon is boson (diquark) in the former. This implies that the ground 
state of the 2-color system at finite baryon density in the color-confined phase is an 
interacting boson system, i.e. a Bose liquid, although the quark Fermi liquid may 
be realized at high baryon density in the color-deconfined phase. How this Bose 
liquid changes its character as a function of the temperature T, the quark chemical 
potential \x and the quark mass m is an interesting question by itself and may also 
give us a hint to understand physics of the color superconductivity in 3-color QCD 
in which the crossover from the Bose-Einstein condensate of tightly bound quark 
pairs to the BCS type condensate of loosely bound Cooper pairs may take placeP 

In this article, we revisit the thermodynamics of the strong coupling limit of 2- 
color lattice QCD with chiral and diquark condensates which was originally studied in 
Ref. HJ). Our main purpose is to analyse its phase structure and the interplay among 
the chiral condensate (qq), the diquark condensate (qq) and the baryon density (q'q) 
as functions of T, \i and m. This would give us a physical insight and useful guide 
for the actively pursued lattice QCD simulations of the same system in the weak 
couplingP 

§2. Formulation 



We derive an effective free energy for meson and diquark fields starting from 
the lattice action with the single component staggered fermion. In the strong cou- 
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pling limit g —* oo, the gluonic part of the action vanishes because it is inversely 
proportional to g 2 . Then the action on the lattice is given by only fermionic part; 

1 d 

S = m^xWxW + g ^2^2vj( x ) {x{x)Uj{x)x{x + j) - x(% + j)uj(x)x(x)j 

X X j = l 

+ 2 {x(x)e ll Uo(x)x(x + 6) - x(x + b)U ] (x)e~" X (x)} , (2-1) 

X 

with 

r!o(x) = 1 , ^(x) = (-l)SLi^-i . ( 2 . 2 ) 

X stands for the quark field in the fundamental representation of the color SU(2) 
group and is the SU(2) valued gauge link variable, d represents the number of 
spatial directions. 

This action has a global Uy(l) x Ua(1) symmetry at /i ^ in the chiral limit 
m = 0. Uv(l) corresponds to the baryon number conservation, which is sponta- 
neously broken by the diquark condensation. On the other hand, the action has a 
larger symmetry U(2) for m = jj, = due to the Pauli-Giirsey's fermion— anti-fermion 
symmetry. This is the special feature of 2-color QCD. 

Here we summarise how to derive the effective free energy for the meson and 
diquark system from the original action Eq. (|2-lj) . 

1. Large dimensional (1/d) expansion is employed in the spatial directions in order 
to facilitate the integration over the spatial link variable Uj. 

2. Bosonization is performed by introducing the auxiliary fields a for XX an d A 
for XX- Then the mean field approximation is adopted for the auxiliary fields. 

3. Integration with respect to x-, X an d t^o are accomplished exactly. 

After performing these procedure,^ we can write down the effective free energy 
for the chiral condensate (a) and the diquark condensate (A) as follows; 

F cS [a,A] = ^ 2 + ^|Z\| 2 -riog{l + 4cosh( J E + /T)-cosh(S_/r)} , (2-3) 
where E± is the excitation energy of (anti) quasi-quarks, 

E± = arccosh (yj (1 + M 2 ) cosh 2 fi + (d/2) 2 \A\ 2 ± M sinh (j\ (2-4) 

with dynamical quark mass M = m + ad/2. In the chiral limit m = with zero 
chemical potential fj, = 0, this is a function only in terms of a 2 + \A\ 2 . As a result, 
the effective free energy is invariant under the transformation mixing the chiral 
condensate with the diquark condensate. This symmetry comes from the Pauli- 
Giirsey symmetry of the original action at m = [i = 0. The computational details to 
derive the effective free energy and its analytical properties can be found in Ref . QJ . 

§3. Numerical Results and Discussions 

We determine the chiral condensate a and the diquark condensate A numerically 
by minimising the effective free energy in Eq. (|2-3|) , The quark density p = —dF e ^/dp 
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Fig. 1. Chiral condensate a, diquark condensate A and quark density p for m = 0.02 with d = 3. 
In the left panel they are plotted as functions of the chemical potential fi at zero temperature. 
In the right panel they are plotted as functions of temperature T for a typical value of \i. All 
the dimensionful quantities are in unit of the lattice spacing, which is implicitly understood in 
other figures. The dotted line indicates a total magnitude of the condensates \J a 2 + A 2 . 



is also calculated. The results are shown in Figs. [I\for small quark mass m = 0.02. 

First we consider the chiral and diquark condensates as functions of chemical 
potential (the left panel of Fig. There exist two critical chemical potentials, 
the lower one p l ° w and the upper one pc P - We can understand the behavior of the 
condensates as the manifestation of two different mechanisms: One is a continuous 
"rotation" from the chiral condensation to the diquark condensation above p = p)? w 
with the total condensate y/ a 2 + A 2 varying smoothly. The other is the "saturation 
effect"; quark density forces the diquark condensate to decrease and disappear for 
large p. 

The "rotation" can be understood as follows: As we have discussed, the effective 
free energy at vanishing m and p has a symmetry between the chiral condensate 
and the diquark condensate. The effect of m (p) is to break this symmetry in 
the direction of the chiral (diquark) condensation favored. Therefore at finite m, a 
relatively large chiral condensate predominantly appears for small p region. Once 
p exceeds a threshold value //j. ow , the chiral condensate decreases while the diquark 
condensate increases, because the effect of p surpasses that of m. 

As p becomes larger, the diquark condensate begins to decrease in turn by 
the effect of the "saturation" and eventually disappears when p exceeds the upper 
critical value pc P (order of unity for T = 0). On the other hand, the quark density 
p increases until the saturation point (p = 2) where quarks occupy the maximally 
allowed configurations by the Fermi statistics. Those behaviors of A and p are also 
observed in the recent Monte-Carlo simulations of 2-color QCDp^ 

Next we consider the chiral and diquark condensates as functions of T (the right 
panel of Fig. At low T, both the chiral and diquark condensates have finite 
values for p = 0.2. The diquark condensate decreases monotonously as T increases 
and shows a second order transition. On the other hand, the chiral condensate 
increases as the diquark condensate decreases so that the total condensate \J o 2 + A 2 
is a smoothly varying function of T. The understanding based on the approximate 
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Current Quark Mass m Fi S- 3 - phase structure of strong coupling 2- 

color QCD in the T-p,-m space. The surface 
Fig. 2. Phase diagram of strong coupling 2- separates the region where A = (outside) 

color QCD in the fi-m plane at T = 0. from the region where A ^ (inside). 

symmetry between chiral and diquark condensates is thus valid. An interesting 
observation is that the chiral condensate has a cusp shape associated with the phase 
transition of the diquark condensate. 

Now we show the phase diagram of the strong coupling 2-color QCD in the 
fj,-m plane at T = in Fig. [21 The lower right of the figure corresponds to the 
vacuum with no baryon number present, p = 0. On the other hand, the upper left 
of the figure corresponds to the saturated system, p = 2, in which every lattice site 
is occupied by two quarks. There is a diquark superfluid phase with A ^ and 
< p < 2 bounded by the above two limiting cases. It is worth mentioning here that 
we can see the corresponding system in the context of condensed matter physics; 
the hardcore boson Hubbard model has a similar phase diagram in which superfluid 
phase is sandwiched by Mott-insulating phases with zero or full density!^ 

Finally, the phase structure in the three dimensional T-p-m space is shown in 
Fig. El The diquark condensate has a none- vanishing value inside the critical surface 
and the phase transition is of second order everywhere on this critical surface. 
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